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EXECUTIVE  SUMMARY 


This  report  describes  a  computer  simulation  of  human  temperature  regulation  for 
predicting  important  physiological  responses  during  exercise  and  exposure  to  heat  stress. 
Although  the  constnrct  is  similar  to  existing  thermoregulatory  simulations,  several  novel 
features  have  been  formulctted.  First,  the  skin  is  not  modeled  as  a  single  compartment 
but  as  two  compartments  corresponding  to  the  deep  vascular  and  the  superficial 
avascular  layers.  Second,  new  algorithms  are  used  to  control  regional  blood  flow  to  skin, 
to  muscle  and  to  core  (wscera  and  other  internal  organs).  Third,  a  new  algorithm  that 
modulates  conductance  between  the  vascular  and  avascular  skin  layers.  Fourth,  there  is 
a  new  model  of  the  central  circulation  that  controls  cardiac  stroke  volume  and  restrains 
blood  flow  to  muscle  and  to  skin  when  demands  exceed  the  maximal  cardiac  output. 
Finally,  the  complexity  of  the  passive  system  is  tailored  to  match  the  problem  of  exercise 
in  hot  environments  where  internal  and  surface  body  temperature  gradients  are  relatively 
small,  symmetrical  and  change  slowly.  Using  a  database  from  7  independent  sturfies  of 
widely  varying  workload  and  environment,  the  simulation  was  validated  by  comparing  the 
standard  deviations  (sd)  of  variable  means  (internal  temperature,  skin  temperature  and 
heart  rate)  in  the  data  sets  with  the  root  mean  squared  deviations  (rmsd)  between  the 
variable  means  and  the  corresponding  simulator  outputs.  Except  for  skin  temperature,  for 
which  there  was  some  uncertainty  about  experimental  conditions,  the  average  sd  and  the 
average  rmsd  were  similar  in  value;  a  difference  of  ±0.03°C  for  internal  temperature,  0 
beats»min  ’  for  heart  rate  and  ±0.78°C  for  skin  temperature.  This  variability  suggests  that, 
on  the  average  and  within  the  range  of  data  tested,  the  present  simulation  is  able  to 
forecast  results  with  the  same  order  of  precision  as  laboratory  studies  of  the  same 
problems.  Possible  military  applications  include  the  accurate  prediction  of  heat  strain  and 
adjustment  of  wotk:rest  cycles  during  intermittent  work  in  protective  clothing  and  in 
confined  spaces  zind  determination  of  optimal  recovery  periods  following  exercise  and/or 
exposure  to  hot  environments.  A  further  application  of  the  simulation  could  be  to  assist 
in  the  design  of  militarily  relevant  studies  on  the  competing  demands  of  work,  heat,  sleep 
deprivation  and  chemical  protective  agents. 

Keywords 

mEtthematical  models;  human  temperature  regulation;  regional  blood  flow;  cardiac  stroke 
volume  control;  skin  thermal  conductance 
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Introduction 

A  computer  simulcition  can  serve  as  a  vector  for  translating  specific  physiological 
knowledge  into  practical  guidelines  for  solving  real-world  problems.  Of  particular  military 
and  industrial  relevance  would  be  a  simulation  of  thermoregulation  designed  specifically 
to  forecast  temperature  changes  within  the  body,  as  well  as  cardiovascular  and  other 
physiological  responses  during  work  in  the  heat.  Such  a  simulation  would  be  a  useful 
tactical  decision  aid  for  predicting  the  level  of  heat  strain  and  casualty  rates  likely  to  arise 
from  different  scenarios  of  activity,  environment  and  protective  clothing  and  for  optimizing 
personnel  turnover  rates  and  adjusting  the  time  periods  in  workrrest  cycles. 

Computer  simulations  of  thermoregulation  typically  incorporate  sets  of  equations 
describing  systems  that  are  physically  passive  and  physically  active  in  heat  transfer 
processes.  The  passive  system  defines  the  geometry  of  the  body  container  space  as  one 
or  more  segments,  divides  the  segments  into  one  or  more  tissue  compartments,  and 
calculates  temperature  distributions  and  rates  of  heat  transfer  w.thin  compartments, 
segments  and  between  the  container  and  the  environment.  The  active  system  defines  the 
response  of  physiological  control  mechanisms  that  attempt  to  alter  rates  of  heat  transfer 
in  the  passive  system  in  response  to  deviations  in  compartment  temperatures  from  certain 
threshold  levels  (31).  Passive  and  active  systems  interact  in  a  closed-loop,  proportional- 
type  control  system  (16). 

The  purpose  of  this  report  is  to  describe  an  uncomplicated  mathematical  simulation  of 
thermoregulation  with  passive  and  active  systems  specificaliy  designed  for  predicting  the 
time-course  of  critical  physiological  responses  of  exercising  humans  during  heat  stress. 
The  following  features  of  this  simulation  are  unique;  (a)  the  skin  is  modeled  as  two 
components:  a  vascular  layer,  where  heat  is  exchanged  with  the  blood,  and  a  superfidal 
avascular  layer,  where  heat  is  exchanged  with  the  environment;  (b)  new  algorithms  are 
incorporated  for  controlling  blood  flow  to  skin,  to  muscle  and  to  core  (viscera);  (c) 
increased  heat  transfer  efficiency  accompanying  cutaneous  venodilation  is  modeled  with 
an  algorithm  for  controlling  thermal  conductance  between  the  vascular  and  avascular  skin 
layers;  and  (d)  a  new  paradigm  of  the  central  circulation  is  proposed  incorporating 
algorithms  for  controlling  cardiac  stroke  volume,  and  for  reducing  regional  blood  flow 
when  demands  for  total  blood  flow  exceed  maximal  cardiac  output.  Some  non- 
thermoregulatory  influences  of  exercise  on  circulatory  components  of  the  active  system 
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are  also^cbn^red  (18^).  The  a'mutation  is  fu^er  evaluated  in  this  repoit=t^ 
comparing  simulator  outputs  with  data  obtained  from  seven  independent  studies  of 
exerose  in  disparate  emrironments. 


Methods 


DEVELOPMENT  OP  THE  SIMULATION 

Simulations  of  thermoregulation  are  usually  constructed  with  spedfic  intentions  and 
definite  applications  in  mind.  The  structural  complexitv  that  must  be  bmit  into  the  passive 
and  active  components  of  a  simulation  largely  depends  upon  the  problem  under 
investigation  (8,10,14,21,31,32).  A  successful  simulation  to  predict  body  temperature 
responses  in  and  around  the  thermal  comfort  zone  has  been  devised  as  one-dimensional 
heat  flow  between  only  two  physical  subdhridons  or  compartments;  a  core  and  a  shell 
(10).  At  the  other  extreme  are  simulations  of  cold  vrater  immersion  or  of  asymmetrical 
heating  and  cooling  during  space  flight  which  reqttire  heat  flow  in  more  than  one 
dimension  between  highly  subdivided  passive  systems  because  of  large  and  rapidly 
changing  temperature  'iradients  across  the  tissues  that  may  not  be  symmetric  throughout 
the  body  (14,  21).  The  optimal  passive  system  of  a  simulation  to  mimic  work  in  hot 
environments  lies  somewhere  in  between  these  two.  On  the  one  hand,  a  human  heat 
strain  simulation  can  employ  fairly  simple  mathematical  forms  without  sacrificing  much 
accuracy  because,  during  exercise  and  environmental  heat  stress,  tissue  temperature 
gradients  are  relatively  small,  can  be  considered  symmetrical,  and  change  rather  slowly 
(14, 32).  On  the  other  hand,  for  accurate  and  comprehenwe  predictioti  of  responses  to 
exerdse  under  heat  stress,  a  simulation  must  model  the  thermal  and  rarculatory  changes 
within  several  body  compartments.  This  requires  a  rather  complex  active  system  and  a 
passive  system  that  is  more  highly  strudured  than  the  core-shell  model. 


The  Passive  or  Controlled  System 
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The  human  body  is  modeled  as  a  single  segment:  an  upright  cylinder  consisting  of  five 
concentric  annular  compartments  (core,  muscle,  fat,  vascular  skin,  avascular  skin)  and 
an  interconnecting  central  blood  compartment  (Rgure  1).  A  one  segment  model  of  the 
passive  systein  was  chosen  overa  more  anatomical  representation  because  data  on  heat 
exchange  coefficients,  clothing  insulation  and  water  vapor  permeability  constants,  and  the 
bulk  of  temperature  and  circulatory  data  needed  to  validate  the  simulation  regard  the  body 
as  a  single  entity,  not  as  muKiple  segments  each  with  its  own  physiological  effector 
responses  and  internal  temperature  distribution. 


210  kg  32.2  kg  9-8  kg  0.7  kg  3.5  kg 


Figure  1.  Cross-section  of  cylindrical  nodel  containing  5  concentric  annular 
tissue  coisparttrents.  Dii&ensions  are  for  an  individual  with  W  «  70  kg,  ”  1'6 
rn*. 

Within  this  upright  cylinder,  heat  is  conducted  radially  between  adjacent  annular 
compartments.  There  is  no  axial  heat  conduction  and  the  cylinder  ends  are  perfectly 
insulated.  Each  compartment  is  homogeneous  and  isotropic  and,  at  any  instant,  of 
uniform  temperature  throughout.  Intercompartmental  resistance  to  heat  flow  occurs  at  a 
single  imaginary  concentric  annular  surface  located  midway  between  the  centers  of  mass 
of  adjacent  compartments.  Between  the  vascular  and  avascular  skin  layers  conductance 
is  controlled  by  the  active  system.  Values  for  thermal  conductance  between  other 
compartments  are  constant  and  vary  only  with  the  anthropometric  characteristics  of  the 
particular  subject. 
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Cylinder  and  compartment  dimensions  are  computed  from  three  anthropometric 
characteristics:  body  weight  (W).  stature  (H)  and  percent  of  body  weight  that  is  fat  (%fat). 
Overall  body  density  (p),  in  g-cm*^.  is  estimated  from  the  empirical  relationship 
p  =  (457)/(%fat-<412),  as  proposed  by  Brozek  et  al.  (5).  This  formula  assumes  that  the 
body  is  made  up  of  only  two  components:  fat  (p^  =  0.9)  and  nonfat  (p^  =  1.1).  V  is 
calculated  as  W/p  and  the  surface  area  (A„),  in  cm^,  is  given  by  the  empirical  DuBois 
formula  (6):  Ao=  (2020)(W)‘''“®(H/100)‘'”®.  The  length  ofthe  cylinder  (L)  is  obtained  from 
the  expression  (Ao)*/(4nV).  and  the  radius  of  the  cylinder  (r)  from  {2V)/(Ao).  The  volume 
of  the  fat  compartment  (V,J  is  obtained  from  (W)(%fat)(pw)-  The  total  volume  of  all  non¬ 
fat  compartments  (V^)  is  calculated  as:  (W)(1C0  -  %fal)(prtJ/{100).  is  then 
distributed  in  the  following  way:  37%  to  com,  54%  to  muscle,  1%  to  vascular  skin,  6%  to 
skin  and  2%  to  the  central  blood  compartment  (30). 

The  procedure  for  calculating  compartment  dimensions  and  themral  conductances 
within  the  passive  system  is  similar  to  that  described  by  Stolwijk  and  Hardy  (30).  Starting 
with  the  outermost  avascular  skin  layer  and  =  V,  the  inner  radius  of  each  annular 
compartment  is  calculated  in  succession  from  the  formula: 


[cm] 
tc'L 

The  radii  of  concentric  cylindrical  surfaces  located  at  the  center  of  mass  of  each 
compartment  am  then  obtained  from: 


(cm)  (2) 

k-L 

The  radius  and  area  of  an  imaginary  concentric  cylindrical  surface  lying  at  the  midpoint 
between  the  centers  of  mass  of  adjacent  compartments  a  and  b  are  obtained  from 
Equations  3  and  4,  respectively: 

V-  =  [cm]  and  V*  = 
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and  length  of  the  conduction  path,  I,  by: 

1^  =  *^r/“*-*  [on]  (5) 

Finally,  intercompartmental  conductance,  1^  is  calculated  from: 


[IV-'C-’J  (6) 

** 

where  K,  and  1^  are  values  of  the  thermal  conductivity  for  compartments  a  and  b  in 
WHaiK^C’^cm^  (f4,30). 


Table  1  Assumed  distributioii  of  body  weight,  mating  heat  production  and  resting  cardiac  ou^t  for  a 
hypothetical  subject  (W  =  70  kg.  A^  s  1 8  m’.  =  141  and  estimated  lumped  values  rf  qredlic  heat,  specific 

thermal  conductivity  and  conductance  for  a  6-nade  cylioderical  represenution  of  this  saliject. 


j  Distribution  of  Body  j 


Fstimafed 


Weight,  Hetalnlism 


EstimaCed 


Thenoal  |  Intcrdlonqiartaient 


Core  ~| 
* 

Muscle 

i 

Fat 

32 

!  “  i 

— — I 
92  1 

51  1 

5.43-10*  j 

1 

30.3 

46 

12  1 

1 

4 

63  ! 

4.18-10-’  1 

23J5 

14 

'  4  : 

1  > 

j  ! 

1  I 

J 

i  42  i 

159-10*  ! 

150.5 

Vascular  9un 

‘  1 

f 

!  0  i 
:  1 

3 

63  i 

4.18-10*  j 

20-400 

AvasetUer 

Skin 

1  ® 

i  >  i 

i  t 
1 

0 

63  j 

2.09-10*  1 

1 

1 

Central  Blood 

1  ^ 

1 

1  0 

0 

i  i 

i 

4.18-10*  } 

A  schematic  dagram  of  the  passive  system  with  computed  dimensions  for  a 
hypothetical  individual  (W  =  70  kg.  A,  =  1.8  m^  is  sfiown  in  Figure  1  and  Table  1  lists 
values  of  passive  system  constants  for  an  individua]  of  this  size. 
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Haat  Prnrilirtinn  Rates  of  metabolic  heat  production  in  the  core,  fat  and  avascular 
skin  compartments  are  assumed  to  be  fixed  percentages  of  total  resting  metabolism  (AO. 
as  shown  in  Table  1  (30).  Heat  production  of  the  muscle  layer  {HJi  is  variable, 
depending  on  total  energy  expenditure  (MJ  and  external  work  performed  (W).  It  is 
assumed  that  the  blood  and  vascular  skin  compartments  produce  no  heat  at  all. 

Heat  Transfer  With  Blood.  All  compartments  (except  the  avascular  skin  layer) 
exchange  heat  with  blood  flowing  through  them  by  convection.  Blood  enters  each  tissue 
compartment  at  the  temperature  of  the  central  blood  compartment  CT^)  and  returns  to  the 
central  blood  compartment  at  the  current  temperature  of  the  tissue  compartment.  Thermal 
mixing  in  the  central  blood  compartment  is  instantaneous.  There  is  no  heat  exchange  virith 
blood  while  in  transit  to  or  from  the  destination  compartment  and  there  is  no  counter- 
current  heat  exchange  between  arterial  and  venous  elements.  Rates  of  blood  flowthrough 
core  (BF„).  through  muscle  (BF.J  and  through  vascular  skin  (BF.J  are  independently 
controlled  by  the  active  system  in  response  to  the  physiological  demands  of  exercise  and 
thennoregulation.  The  rate  of  blood  flow  through  fat  (BFJ  is  considered  to  be  small 
(1.2%  of  basal  cardiac  output)  and  unchanging. 

Heat  Transfer  With  The  Environment  Heat  exchange  between  skin,  clothing  and 
environment  is  assessed  using  the  partitional  calorimetric  model  of  Gagge  and  Nishi  (1 1 ). 
Heat  is  lost  from  the  skin  surface  by  the  combined  sensible  routes  of  radiation  and 
convection  (R+C)  and  by  the  insensible  route  of  sweat  evaporation  (£).  Heat  loss  also 
occurs  from  the  cemral  blood  compartmem  via  the  respiratory  passages  through 
convection  and  evaporation  (C^+E^. 

Sensible  Heal  Loss.  (R+C)  is  given  by: 

(R*C)  =  /v,h(VTJ.  [W-m^  ^ 

where  h  is  the  combined  coefficient  for  heat  transfer  by  radiation  and  convection,  defined 
as  h  =  (h,  +  hj.  T„  is  operative  temperature,  defined  as  T,  =  (h,T,+  h  JJ/(h).  and  F'^  is 
a  dimensionless  climatic  efficiency  factor,  defined  as  1/(1  +  0.155l*h^.  I^o  is  the  clothing 
insulation  value  in  Burton’s  do  units.  T,  is  the  mean  rarSant  temperature  (“C)  and  T.  is 
ambient  air  temperature  (“C). 
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h,  is  obtained  from: 


h,  =  Aa(A/A^f, 


flEfUR. 
“t  2 


273 


where 


(8) 


a  is  the  Stefan-Boltzman  Constant  5.67x10**  [W- 
[A/Af)  is  the  fraction  of  body  area  exposed  to  radiation  (C.72),  and 
f^  is  the  Breckenridge  clothing  area  factor,  defined  as  (1  +  0.1  sy  [nd]. 


T^,  the  temperature  of  the  outer  clothing  surface,  is  estimated  by  the  expression: 

s  T,  -I-  -  TJ.  Formuias  used  by  the  program  to  caiculate  the  convective  heat 

transfer  coefficient,  h^  are  shown  in  Table  2.  The  variables  v^  and  are  speeds  of 
air  around  a  stationary  body  and  of  a  body  through  stationary  air,  res|>ectively.  in  m*sec''. 


Table  2.  Formolaa  for  cueaUtmg  h„  the  oomreetiee  beet  earhange  coefficient.  Adapted  fiom  Gagge 


andNishi(ll). 

ActfvftTLnrrd 

FormcdafarA. 

Sittiiig  Rot,  Statranaij  Ann  Work  . 

h,=  lL6(»JI“ 

Sutionaijr  Ticadmill  Walking 

h.  =  6.5(T«»y"  + 

Free  Walking 

h,  =  ♦  lS6(v*T* 

StatkmazT  Bicjrde  at  50  rerqnin'’ 

h.=  5.5  + 

Insensible  Heat  Loss.  E  is  determined  by  the  rate  of  sweat  secretion  (rii„)  and  the 
maximal  rate  of  evaporative  heat  loss  from  a  fully  wetted  skin  surface  is  a 

function  of  the  vapor  pressure  gradient  between  the  fully  wetted  skin  surface  and  the  air 
(P^  -  PJ,  the  evaporative  heat  transfer  coefficient  (hj  and  i,^  Woodcock's  dimensionless 
feictor  for  permeabiEty  of  water  vapor  through  clothing.  The  evaporative  heat  transfer 
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coefficient,  h.,  is  directly  related  to  the  convective  heat  transfer  coefficient,  he  by  the 
Lewis  Relation  (13). 


When  evaporation  is  not  restricted  by  clothing  or  the  environment  then  ((h..X)/Ao) 

in  W«m^,  where  is  in  g>min''  and  X  is  the  heat  of  vaporization  for  sweat  at  35°C: 
40.8  W>min*g'\  The  expression  for  under  conditions  where  evaporation  of  sweat  is 
restricted  is  (13): 

=  forE^&m^  (9) 

where  is  the  modified  dimensionless  [no] 

Nishi  permeation  effidency  factor  =  ^  f^iin,  is  the  Dubois 
surface  area  and  /i,  =  2.2/»^  (af  sea  level)  {W-m'^-Torr'\. 


P,^  is  related  to  by  the  Antoine  Equation  (13): 


=  exp  18.6686 


4030.183 


235 


[Toni 


(10) 


Respiratory  Heat  Loss.  {C„+E^  is  directly  related  to  ventilation  rate  which,  in  turn, 
is  directly  related  to  aerobic  exercise  intensity  {M^  up  to  maximal  levels.  The  combined 
equation  for  convective  and  evaporative  respiratory  loss  is  taken  from  Fanger  (9): 

=  M„[0.0014(34  -  7J  ♦  0.0023(44  -  PJ]  \Wm^  (11) 
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Compartmental  Heat  Balance.  The  usual  conservation  of  energy  statements  apply 
and  have  been  described  in  detail  by  others  (28,  32).  The  general  differential  equation 
describing  the  rate  of  change  in  heat  content  of  the  n**  annular  compartment  at  any  time 
t  is: 

=  H„(t)  *  iK^aTat)-Tjn)n  - 

The  spedfic  equations  for  each  solid  compartment  are: 

4°-  =  lK.JiT^-TJ]  -  lBF^P^-cJiT„-TJi  (12a) 

*  IK^JiTc-TJ]  -  IK^T^-TJ]  -  IBF^-p^-cJiT^-TJ]  (12b) 

♦  IK^T^-TJ\  -  IK^JiT^-T^]  -  Vrj]  (12r^ 

^  -  [BF«»-P«cjr^-Tj]  (12d) 

-  AJi(R*C)  *  q  (12e) 

and  for  the  central  blood  compartment: 

=  P-c^iKVrjBFj  ♦  l(r„-rjBFj  ♦  kvtjbfj  ♦ 

+  [(T-^-TJBFJ)  -  iC^^EJ  (120 

Each  cfifferential  equation  is  evaluated  by  a  simple  numerical  integration  procedure  that 
is  based  on  the  assumption  that  if  at  is  made  very  small,  then  the  value  of  dQ^dt  at  time 
t.  is  appron'mately  the  same  as  the  value  of  ACVAt  over  the  interval:  (V^)  to  U  0  0,31 ). 
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Then  the  change  in  heat  content  of  compartment  n,  AQ„,  is  obtained  ua'ng  the  following 
approximation: 

AQ„(UlU  =  |Q«{U  Af  [IV-minl  (13) 

From  this,  AT„(tJ  is  calculated  from: 

ATJ/tJ  =  [C»  or  K[ 

"  P.c.V-. 

(14) 

where  p„,  c„  and  V„  are  the  density,  heat 
capacity,  and  volume  of  the  rtth  compartment 

At  any  chosen  time  of  interest  (t,)  the  temperature  of  the  compartment  is  obtained  by 
iteration: 


Ut)  =  \  *  EaW,  rq 
*>1 

where  i  =  —  +  ^ 

Af 

and  T.  is  the  initial  temperature  condition 


There  would  rarely  be  reason  to  choose  successive  values  of  t|  less  than  0.25  min 
apart  and  At  is  initially  set  by  the  program  at  0.025  min,  thus  providing  a  minimum  of  10 
iterations  per  reporting  interval.  Under  most  circumstances  this  iteration  rate  is  adequate 
to  satisfy  the  assumptions  regarding  the  use  of  Equation  13.  However,  if  AT„  >  0.1  C° 
during  any  interval  then  At  is  automatically  reduced  by  the  program  for  the  remainder  of 
interval  t,  according  to  the  formula: 


(0.1){AU 


At  is  automatically  reset  to  0.025  min  for  the  beginning  iteration  in  the  next  time  of  interest 
period  (tj. 
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Active  System 


The  active  system  consists  of  mathematical  algorithms  that  describe  the  behavior  of 
physiological  effectors  whose  outputs  attempt  to  control  internal  temperature  by 
redistributing  blood  flow  and  among  tissue  elements  and  by  increasing  the  potential  for 
environmental  heat  loss  by  increasing  sweat  secretion.  On  occasion,  increased 
thermogenesis  from  shivering  may  be  activated.  Inputs  to  the  simulation's  active  system 
include  temperature  levels  within  the  passive  system,  primarily  from  the  central  blood 
compartment  (T^)  and  from  skin  (t^  and  also  inputs  from  age,  energy  expend'rtu.ro  (M„ 
or  VoJ  and  from  heat  production  (H^J.  Effectors  include  rates  of  blood  flow  and 
convective  heat  flow  among  four  of  the  five  tissue  compartments,  conductance  between 
vascular  and  avascular  skin  compartments,  the  rate  of  sweating  and  the  rate  of  shivering. 
Also  the  level  of  cardiac  stroke  volume,  whidi  is  influenced  by  both  work  intensity  and 
thermal  factors,  is  controlled.  Active  system  algorithms  and  values  of  heat  exchange 
between  skin,  clothing  and  the  environment  are  recalculated  during  each  iteratio.i. 


Chanae-ot-State  Lags 


Discontinuities  in  physiological  responses  occum'ng  during  the  transition  between 
exercise  levels  can  be  particularly  troublesome  when  there  are  large  changes  in  workload. 
First-order  lags  are  introduced  so  as  to  restrain  large,  sudden  changes  in  calculated  blood 
flows  and  stroke  volume  that  would  otherwise  accompany  change-of-state  discontinuities. 
These  lags  are  described  by: 


-  XJ 


( 


1  -  exp 


-0.693 L 


P 


where 


(17) 


X„.,  is  the  new  time-lagged  value  for  variable  X,  X,„  is  the  time-lagged  value  for  X  during 
the  last  iteration  intenral,  is  the  new  non-lagged  value  for  X  as  calculated  by  the 
algorithm,  t„  is  the  elapsed  time  into  the  current  time  of  interest  period  and  tg^  is  the 
response  half-time  (from  30  sec  to  2  min,  depenrSng  upon  the  function). 
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Skin  Blood  Flow 


The  rate  of  thermoregulatoiy  skin  blood  flow  (BF^  is  described  primarily  as  a  linear 
function  of  change  in  hypothalamic  temperature  from  a  reference  point,  but  is  also 
modulated  by  skin  temperature  (t^.  by  posture,  by  exercise  intensity  and,  transiently,  to 
the  initiation  and  cessation  of  exercise  (4, 17-19, 23).  The  algorithm  used  in  the  present 
simulation  for  control  of  BF.^  is  based  primarily  on  an  empirical  regression  model 
suggested  by  Roberts  and  Wenger  during  upright  exercise  and  heat  stress,  but  also 
incorporates  effects  of  activity  level  as  inferred  from  other  work  (18, 23). 


Figure  1,  Algorichjs  for  calculating  Haxiisua  is  determined  by  activity 

(e),^  threshold  by  activity  tt  t„  and  level  by  Lags  introduced  by  Eq.  17. 


A  quantitative  representation  of  this  control  algorithm  is  rfiagrammed  in  Rgure  2.  BF,^ 
is  modeled  as  a  linear  function  of  central  blood  temperature  (Tu).  Both  ^nd  exerdse 
shift  this  function's  intercept  on  the  Ty  axis:  increasing  reduces  the  value  of  T^ 
necessary  to  produce  a  given  value  of  BF^  while  exercise  increases  the  value  of  T^ 
required  to  produce  a  given  value  of  BF.^  Exercise  also  produces  graded  reductions  in 
maximal  skin  blood  flow:  from  7.0  Fmin  ’  at  rest  to  5.0  l•min■’  at  a  Vo,  of  2.0  Umin  ’  (18). 
Minimum  BF„,  during  intense  vasoconstriction  is  set  at  30  ml«min'\  The  assumption  is 
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made  that  all  'rest'  is  in  the  upright  seated  position  and  that  all  “woik*  is  in  the  upright 
standing  position. 


Vascular-Avascular  Skin  Conductance 

Cutaneous  veins  ate  quite  distensible  and  fill  with  blood  as  BF^  increases  (26).  An 
increase  in  transit  time  of  blood  through  the  cutaneous  venous  bed  would  presumably 
accompany  this  venous  distension,  resulting  in  greater  heat  transfer  between  blood  and 
skin.  Effectively,  then,  the  value  for  conductance  between  these  layers  {K.,1,^  increases 
as  BF,^  increases.  To  date,  there  is  no 
direct  evidence  to  document  the  nature 
of  the  relationship  between  BF^  and 
The  presumption  is  made  here 
that  each  ten-fold  increase  in  BF,^ 
increases  by  a  fixed  amount 
(19).  The  maximal  value  for 
obtained  rationally  by  calculation,  as¬ 
suming  maximal  vasodilation  and  per¬ 
fect  thermal  contact  (Equation  6  and 
Table  1).  The  minimum  value,  10 
cm*»m'*»min'’,  assumes  maximal  va¬ 
soconstriction,  so  that  for  heal  transfer 
purposes  the  vascular  skin  layer  is 
essentially  bloodless  and  similar  to  the 
avascular  skin  layer.  K,riui>,  modeled 
as  a  function  of  log  BF,„„  is  shown  in 
Rgure  3. 


Core  Blood  How 

While  the  core  contains  the  skeletal  system,  the  nervous  system  and  connective 
tissue,  these  components  make  negfigible  demands  on  carcfiac  output.  The  core 


Figortt  3.  Hypothetical  relationship  between 
skin  blood  flow  and  skin  conductance. 
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components  receiving  the  largest  blood  flows  at  rest  are  kidneys  and  splanchnic  regions 
which  obtain  over  a  third  of  the  resting  cardiac  output  With  the  addition  of  exercise  and/or 
heat  stress,  blood  flow  to  splanchnic  and  renai  vasculatures  can  be  reduced  by  as  much 
as  70%  (26). 


Figure  4«  BF,,  is  calculated  iron  an  observed  empirical  association  with  heart 
rate.  Basal  BFcr  assuned  to  be  92%  of  basal  cardiac  output.  Frc»a  Howell  (25). 


The  paradigm  for  controlling  core  or  visceral  blood  flow  (BFJ  is  based  on  a 
regression  modei  of  splanchnic  and  renal  blood  flow  on  heart  rate  (Rgure  4)  that  stems 
from  many  observations  taken  during  rest  and  exerdse  with  and  without  heat  stress  (25). 
The  data  fall  on  two  nearly  parallel  lines  (one  for  rest  and  another  for  exercise)  regardless 
of  whether  heart  rate  was  elevated  b>  heating  or  by  lower  body  negative  pressure. 
AKhough  it  is  unlikely  that  decrements  in  are  linked  causally  to  increments  in  HR,  the 

operational  assumption  is  made  that  they  are  both  well-related  to  a  third  variable  that  is 
independent  of  bltx>d  flow. 
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Muscle  Blood  now 

This  model  for  the  control  of  BF^  is  a  function  of  (a)  the  oxygen  requirements  of 
skeletal  muscle  (MuVOj),  (b)  the  demands  for  total  cardiac  output  (CO^),  (c)  the  maximal 
cardiac  output  (CO^  and  (d)  the  level  of  oxygen  extraction  from  the  blood  (PclOjext). 
In  general,  BF„„  is  assumed  to  be  a  linear  function  of  MuVo,;  PctO^ext  is  assumed  to  be 
60  percent  and  CO,^  is  assumed  to  be  less  than  CO„„  (Rgure  5).  If  CO,^  at  anytime 
becomes  greater  than  CO^,  then  PctOjext  is  incrementally  increased  to  force 
compensatory  reductions  in  BF^^  The  change  in  PctOjext  is  based,  in  part,  on  a 
hypothetical  scheme  proposed  by  Rowell  (26).  As  shown  in  Rgure  5,  MuVo,  can  be 
impacted  by  both  activity  level  and  shivering  thermogenesis. 


Figure  5.  BF.,  is  a  function  of  cuscle  Oj  consumption.  Percent  extraction 
variable  depending  on  competing  demands  for  cardiac  output.  Shivering  model 
from  (31). 
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Cardiac  Stroke  Volume 


Stroke  volume  (SV)  is  a  complicated  function  whose  value  is  influenced  by  many 
factors.  Only  two  are  considered  here;  the  level  of  exertion  and  the  level  of  skin 
temperature,  as  shown  in  Figure  6.  Starting  at  rest,  assumption  of  upright  posture 
displaces  central  intravascular  volume  to  the  periphery,  thus  reducing  filling  pressure  and 
SV.  At  this  point  maximal  stroke  volume  (SV„^  is  at  its  lowest,  about  85  cm^  (26). 
Beyond  this,  (SV„,„)  increases  oirectly  with  activity  level,  up  to  a  t/0,  of  about  2  ^min■' 
where  it  attains  a  high  of  about  130  cm’  (2). 

When  is  high,  SV  can  be  diminished  because  of  displacement  of  central 


Figure  6.  Uprighc  cardiac  scroke  volune  is  modeled  as  a  function  of  the  level 
of  exertion  and  of  skin  temperature.  Lag  function  described  by  Eg.  17. 


intravascular  volume  into  cutaneous  veins  (25).  This  aspect  of  SV  control  is  modeled  in 
the  following  way.  At  rest  this  displacement  effect  is  nil;  the  value  of  SV  is  equal  to  resting 
SV„a.  regardless  of  the  magnitude  of  T*.  As  energy  expenditure  and  SV„,„  increase,  the 
impact  of  increasing  becomes  more  pronounced  and  the  final  value  of  SV  can  be  as 
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much  as  25  ml  lower  than  when  is  cooL  It  is  cautioned  that  the  quantitative  aspects 
of  this  strategy  are  largely  hypothetical. 


Cardiac  Output  and  Heart  Sate 

At  any  instant  the  actual  cardiac  output  (CO)  is  the  smaller  of  required  cardiac  output 
(CO^  and  maximal  cardiac  output  (CO-,.).  Heart  rate  (HR)  is  a  variable  determined 
simply  by  the  ratio  of  present  values  for  CO  and  SV.  CO,„  is  determined  by  the  sum  of 
required  blood  flows  to  all  compartments.  CO^  is  determined  by  the  product  of  maximal 
heart  rate  (HR—)  and  the  current  value  for  SV  (Rgure  7).  HR^  is  computed  from  the 
approximate  guideline:  HF)^  =  220  -  age  (1). 


Figure  7*  Relationship  between  CO,m  Rate  of  muscle  O3  extraction  can 

be  increased  by  excessive  demand  for  cardiac  output. 
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Sweating  Rate  and  Shivering 

The  sweating  rate  algorithm  Is  that  proposed  by  Nadel  et  aJ.  (22)  as  a  function  of 
hypothalamic  temperature  (Ty)  and 


(18) 


exp 

10 

Threshold  temperatures  T,^  and  for  T^  are  set  at  36.96°C  and  SS-OO’C,  respectively. 


Figure  t.  Diagra^atic  represer.taciop  of  the  sweating  rate  aigorithn  (Eq. 
18).  Fron  Nadel  ec  ai.  (22). 


Although  this  emulation  is  not  designed  for  use  during  obvious  cold  exposure, 
shivering  can  occur  under  room  temperature  conditions  when  persons  are  recovering  from 
heavy  exerdse  or  heat  exposure.  Therefore  a  shivering  algorithm,  originally  proposed  by 
Stolwijk  and  Hardy  and  shown  in  Rgure  5,  was  included  in  the  simulation  (31). 
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Model  Implementation 


A  flow  chart  summarizing  the  basic  elements  of  program  flow  is  shown  in  Rgure  9.  The 
simulation  was  developed  as  an  interactive  program  and  compiled  in  a  structured  BASIC 
dialect  (Turbo  BASIC  1 .1 ,  Borland  International,  Ina.  Scotts  Valley,  CA)  but  any  common 
programming  language  should  be  satisfactory.  Program  inputs  include  subject 
anthropometric  data  (stature,  weight,  age,  percent  body  fat),  clothing  insulation  and  water 
vapor  permeabifity,  environmental  variables  (air  temperature,  humidity  and  movement), 
type  and  intensity  of  exercise  (leg  work,  arm  work  or  rest;  walking  speed  and  grade  or 
energy  expenditure),  and  length  of  time  cycles. 


Validation  of  the  Simulation 


Data  from  seven  studies  emphasizing  various  aspects  of  exercise  under  heat  stress 
were  chosen  to  test  the  simulator.  Selected  studies,  representing  different  investigators 
and  laboratories,  presemed  enough  detail  about  subjects,  measurement  methods,  activity 
levels,  and  environmental  conditions  so  that  the  protocol  could  be  accurately  input  to  the 
simulator.  Six  of  the  studes  included  tabular  data  and  the  standard  deviation  statistic  or 
included  individual  subject  data  so  that  this  statistic  could  be  calculated.  Data  from  none 
of  these  studies  had  been  used  in  development  of  the  simulator. 

Data  Set  1.:  Ekblom  etal.  (7).  Three  well-trained  subjects  ctothed  in  shorts  exerdsed 
by  continuously  cycGng  at  60%  Vo,max  in  a  comfortable  (21  °C)  environment.  Air 
movement  was  rapid  to  insure  free  evaporation.  Data  for  T„,  and  HR  were  available 
at  50  min  of  exerdse  and  HR  and  were  compared  to  simulator  outputs  programmed  for 
that  time.  This  data  set  tests  the  simulator  response  to  exertion  without  external  heat 
stress. 

Data  SetZ:  Gonzalez  etal.  (12).  Frve  subjects  clothed  in  shorts  each  exerdsed  in  16 
experiments  on  a  bicycle  ergometer  at  28%  Vp,max.  In  each  experiment  subjects  were 
exposed  to  a  different  combination  of  temperature  and  humidity  chosen  to  produce  one 
of  six  effective  temperature  (£7)  states.  Tabulated  steady-state  values  lor  T„,  T^and  HR 
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Flow  chart  of  main  program  elements  In  the  simulation.  Calculated  values  are  reported  at  multiples 


at  40  min  were  compared  with  simulator  outputs  at  the  same  time.  This  data  set  tests  the 
simulation  with  respect  to  a  wide  range  of  environments  without  the  complicating  effects 
of  changing  workload. 

Data  Set  3.:  Kraning  and  Gonzalez  (20).  Four  subjects  exerdsed  in  a  warm 
environment  (30°C.  25%  RH)  for  60  min  under  four  protocols;  continuous  treadmill 
walking  (1.34  m-s  ',  3%  grade)  while  wearing  shorts  or  heavy  semi-permeable  clothing 
and  repeated  cycles  of  walking  (1 .34  m-s'\  0%  grade  for  4  min),  jogging  (2.24  m-s'\  5% 
grade  for  2  min)  and  rest  (seated  for  4  min)  while  again  wearing  either  shorts  or  heavy 
semi-permeable  clottung.  This  data  set  provided  measurements  to  illustrate  the  time 
course  of  responses  to  mild,  fully  compensable  external  heat  stress  and  to  severe, 
uncompensable  external  heat  stress  during  two  different  patterns  of  internal  heat 
production.  Data  were  reported  at  15  second  intervals  for  T„,  and  for  HR  and 
averaged  across  subjects.  Each  averaged  15  second  data  point  was  compared  with  its 
corresponding  simulator  output  over  a  60  min  period. 

Data  Set  4.:  Pandolf  et  at.  (unpublished  obsenrations).  10  subjeds  wearing  shorts 
rested  for  10  min  in  a  hot,  dry  (49°C,  20%  RH)  environment  then  worked  in  the  same 
environrrrent  on  a  stationary  ergometer  at  a  workload  requiring  53%  Vo^ax  for  50  min 
more.  This  data  set  illustrates  the  physiological  consequences  of  heavy  work  under 
desert-Ske  condWons  and  is  typical  of  responses  observed  on  the  first  day  of  heat 
acclimation.  As  in  Data  Set  3.,  T.  and  were  compared  with  simulator  outputs  every 
15  seconds  for  60  min.  HR  data  were  available  for  comparison  at  10  min  intervals. 

Data  Sets.:  Rowell  et  at.  (24).  Average  central  drculatory  and  thermal  responses  of 
6  sui^ects  exercising  at  31%  Vo,max  and  5  subjects  exerddng  at  54%  Vo^ax  were 
measured  while  T^  was  driven  between  28°C  and  40'’C  with  a  temperature  controlled, 
water  perfused  suit  In  simulating  this  data  set  was  made  a  forcing  function  rather  than 
a  dependent  variable,  and  equations  describing  heat  exchange  with  the  environment  and 
clothing  were  eliminated,  was  independentiy  calculated  during  each  iteration  interval 
using  a  set  of  best-fit  polynomial  regression  equations  that  duplicated  the  sidn 
temperature  response  of  the  original  data  set  Using  this  data  set  it  was  possible  to 
compare  directly  data  and  simulation  outputs  for  T„„  and  for  HR,  SV  and  CO  at  5  min 
intervals. 
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Data  Set  6.:  Saltin  ei  at.  (27).  Fourtrained  subjects  exerdsed  in  separate  experiments 
at  either  27%,  46%  or  72%  Vo,max  in  an  ambient  temperature  of  either  10®C,  20®C,  or 
30°C.  Reported  data  for  T„,  t*  and  for  HR  at  the  end  of  each  exposure  (50  min)  were 
compared  with  simulator  outputs. 

Data  Set  7:.  Skokterstrom  (29).  Eght  fire  fighters  each  exercised  by  walking  for  one 
hour  at  2  workloads  (21%  and  31%  Vcynax)  in  4  experiments  while  wearing  either 
(xiveralls  or  fire  fighting  turnout  gear  in  a  15“C  or  a  45“C  environment.  Reported  end 
values  for  T„,  f*  and  for  HR  at  60  min  were  compared  with  simulator  outputs. 


Statistical  Method 


A  summary  statistic,  the  root  mean  squared  deviation  (rmsd),  was  proposed  by  Haslam 
and  Parsons  to  compare  values  of  simulated  variables  with  average  values  of  experimen¬ 
tal  variables  at  corresponding  time  points  (15).  The  rmsd  statistic  is  defined  as: 


rrasd  = 


where: 


d,  =  difference  between  observed  and  predicted  at  each  time  point 
n  =  number  of  comparison  time  points. 


This  statistic  is  comparable  with  the  average  standard  deviation  (sd)  of  the  data  set  The 
ave  sd  wets  available  for  all  but  Study  No.  5. 
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Results 


Comparisons  of  values  for  the  sd  of  experimental  data  and  the  rmsd  between 
experimental  data  averages  and  simulator  outcomes  are  detailed  in  Table  3  and 
summarized  in  Figure  10.  Fbrthree  of  the  studies  (No.  1,  No.  6  and  No.  7)  each  reported 
value  for  the  sd  and  the  rmsd  repre¬ 


sents  a  single  end-point  comparison  at 
each  condition.  In  one  study  (No.  2) 
each  value  of  the  sd  and  the  rmsd 
reported  represents  the  average  of  2 
to  3  end-points  at  each  oondriion.  In 
two  studies  (No.  3  and  No.  4)  each 
reported  value  of  the  sd  and  of  the 
rmsd  represerris  the  average  of  com¬ 
parisons  at  15  second  intervals  for  the 
entire  rime  course:  240  and  200  com¬ 
parisons,  respectively.  In  one  study 
(No.  5)  the  repotted  average  rmsd  is 
based  on  21  comparisons  during  the 
rime  course  of  the  experiment  and 
values  for  the  sd  were  not  available. 


RgurvlO.  Comparison  of  data  standard  deviations  with  tms 
deviations,  sd  not  avaaaUe.  controlled. 
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Internal  Temperature 


In  1  study  (No.  5)  internal  temperature  (!,„)  was  taken  as  blood  temperature  in  the 
right  atrium  (TJ.  in  2  studies  (No.  2  and  No.  4)  as  esophageal  temperature  (T^  and  in 
the  remaining  4  studies  as  rectal  temperature  (TJ.  The  overall  average  sd  for  internal 
temperature  was  0.25°C.  The  rmsd  statistic  for  Tj^  vs.  T„  was  generally  smaller  than  the 
rmsd  statistic  forTj,  vs.  T^  (Rgure  10).  The  average  rmsd  for  T„  vs.  T„  was  0.28'’C  and 
the  average  rmsd  for  T^,  vs.  Ty  was  0.42°C. 


Mean  Skin  Temperature 


Study  No.  5  was  not  included  because  T^  was  controlled  as  an  independent  variable. 
The  average  sd  for  was  0.55°C  while  the  average  rmsd  statistic  was  1.33°C.  In  4  of 
the  6  studies  the  sd  and  rmsd  were  in  close  agreement  but  in  2  studies  (Nc  1  and  No. 
6)  there  was  substantial  cOsagreement;  the  simulation  consistently  overpredicted 
(Rgure  10).  The  subjects  of  these  two  studies  were  reported  to  be  in  a  high  state  of 
physical  training.  Since  physical  training  confers  a  degree  of  heat  acclimation,  perhaps 
these  individuals  had  an  earlier  onset  of  sweating  and  a  lower  T^  than  would  be  expected 
for  non-acclimated  subjects. 


Heart  Rate 


The  value  of  the  average  rmsd  statistic  for  mean  heart  rate  vs  simulator  output  was 
the  same  as  the  value  of  the  average  sd  (12  beats*min'')  although  both  the  sd  and  the 
rmsd  were  quite  variable  between  studies. 
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Sample  Plots  of  Simulator  Outputs 


The  sd  and  rmsd  error  statistics  are  standard  metrics  for  comparing  the  overall 
precision  of  the  simulator  with  the  obsen/ed  variance  of  experimental  data;  however, 
taken  alone  such  comparisons  fail  to  identify  problems  of  tracking  or  consistent 
underprediction  or 
overprediction. 

Correct  tracking  is 
particularly  impor¬ 
tant  when  proto¬ 
cols  include  cha¬ 
nges  in  energy 
expenditure  or 
environmental 
conditions.  Figure 
11  shows  the  time 
course  of  ten 
selected  simulator 
outputs  and  three 
experimental 
counterparts  for 
Study  4.  Panel  a 
shows  predicted 
heat  production  at 
rest  and  while 
cycling.  Panel  b 
compares  the 
time  course  of  T„ 

±  sd  with  simula¬ 
tor  outputs  for  T„ 
and  Tu.  Both 
simulator  temper- 


-10  0  10  30  30  40  $0  00  -10  0  10  20  30  40  50  60 

Time  (minutes) 


Figure  11 .  Time  couise  or  selected  simulator  oulpuls  and  3  enperiniental  counterparts 
for  Study  No  4. 
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atures  overpredicted  the  value  but  was  better  both  in  quantitative  agreement  and 
in  tracking.  Panel  c  compares  the  simulation  output  and  data  for  average  t,  ±  sd. 
Although  the  simulation  tracked  the  features  of  the  data  average,  its  course  was  quite 
exaggerated.  Quantitatively,  however,  the  rmsd  was  slightly  smaller  than  the  data  sd. 
Panel  e  shows  average  HR  data  ±  sd  collected  at  10  min  intervals  compared  with  15 
second  output  from  the  simulator.  Although  the  simulation  consistently  underpredicted  HR, 
time-course  traddng  was  reasonable. 

Figure  12 
shows  the  time 
course  of  the 
same  variables  for 
the  fourth  protocol 
of  Study  3:  a 
highly  stressful 
intermittent  work 
routine.  Panel  a 
shows  the  time 
course  of  predict¬ 
ed  heat  produc¬ 
tion  corresporKfing 
to  cycles  of  walk¬ 
ing,  jogging  and 
recovery.  Panel  b 
compares  the 
time  course  of 
average  T„  ±  sd 
with  simulator 
outputs  for  T„  and 
T^,.  Although  T„ 
slightly  underesti¬ 
mated  T„,  track¬ 
ing  was  good  and 
small  oscillations 
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Figure  12.  Selected  Emulator  outputs  and  3  experimental  counterparts  from  one 
segment  of  Study  No.  3. 
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of  the  two  temperatures  with  workload  were  in  phase  with  each  other  and  similar  in 
magnitude.  The  average  rmsd  was  0.1 1°C  compared  to  an  average  data  sd  of  0.29°C. 
On  the  other  hand,  there  was  poor  correspondence  between  Ty  and  T^.  The  average 
rmsd  was  0.70°C  and  Ty  was  consistently  lower  in  value  than  either  T„  or  T^. 
Qualitatively,  Ty  was  quite  different  than  T,,,  showing  a  three  component  waveform  with 
crests  and  troughs  that  were  more  peaked,  larger  in  magnitude,  and  out  of  phase  with 
Ty.  This  discrepancy  was  not  entirely  unexpected  since  heat  exchange  between  the 
muscle  and  the  central  blood  compartments  is  affected  directly  and  immediately  by  large 
swings  in  muscle  blood  flow  accompanying  changing  exertion  level  (panel  g). 

Panel  ccompares  the  time  course  of  average  ty,±  sd  with  the  corresponding  simulator 
output.  Average  rmsd  (0.56°C)  was  slightly  greater  than  average  data  sd  (0.41  °C). 
Qualitatively,  pulsations  in  data  and  simulator  output  accompanying  changes  in  workload 
were  similar  in  phase,  waveshape  and  magnitude. 

Panel  e  compares  the  time  course  of  average  HR  ±  sd  and  the  simulator  output  for 
HR.  Waveshape  and  phcise  relations  were  similar  for  data  and  simulation.  On  the  whole, 
an  ave  sd  of  10  bpm  compared  with  an  rmsd  of  18  bpm.  Reasonable  magnitude 
agreement  was  achieved  during  heavy  exercise  and  recovery  phases.  However,  during 
the  moderate  exercise  periods  the  simulator  always  •  derpredicted  the  level  of  average 
HR.  Since,  in  this  simulation,  HR  is  not  controlled  but  is  simply  the  quotient  of  CO  and 
SV,  the  error  must  lie  either  in  underestimating  the  level  of  CO  or  in  overestimating  the 
level  of  SV  or  in  both.  Since  neither  CO  or  SV  were  measured  in  this  study,  the  source 
of  error  cannot  be  determined. 

Figure  13  is  a  plot  using  data  from  the  second  protocol  of  Study  5  where  subjects 
exercised  continuously  and  ty,  was  independently  controlled  with  a  water-perfused  suit. 
In  the  simulation,  control  of  ty,  was  achieved  by  dividing  the  skin  temperature  data  curve 
into  4  time  segments,  finding  the  best-fit  5th  order  polynomial  expressions  for  each 
segment  and  then  substituting  these  4  equations  for  the  previous  definition  of  ty, 
(Equation  15).  This  procedure  gave  good  agreement  between  the  experimental  and 
simulated  ty,  forcing  functions  (Panel  c).  Panel  a  shows  that  simulator  heat  production 
was  constant  except  for  a  small  increase  from  shivering  (<  1  %)  when  ty,  fell  below  32°C. 
Panel  b  shows  the  time  course  for  average  Ty,  measured  in  the  right  atrium,  forTy,  and 
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for  simulator  generated  and  Both  simulator  outputs  tracked  average  reasonably 
well  but  both  outputs  overshot  peak  average  observed  during  skin  heating  by  a  large 
amount.  0.5°C.  Overall,  tracking  between  average  and  simulated  was  somewhat 
better  than  between  average  Ty  and  simulated  T„.  T„  was  more  sluggish  than  average 
Ty  or  either  of  the  simulated  temperatures,  particularly  during  rapid  skin  cooling. 

Tracking  between  simulated  and  average  CO  (panel  g}  was  also  reasonable,  especially 
during  skin  heat¬ 
ing,  although 
simulated  CO 
reached  an  equi¬ 
librium  level  earli¬ 
er.  However, 
during  rapid  skin 
cooling,  average 
CO  declined  more 
slowly  than  simu¬ 
lated  CO.  At  the 
same  time,  simu¬ 
lated  SV  rose 
much  faster  and 
to  a  greater  extent 
than  average  SV 
(panel  f).  The 
combined  effect  of 
underpredicting 
CO  and  overpre¬ 
dicting  SV  during 
rapid  skin  cooling 
was  that  simulat¬ 
ed  HR  (Panel  e) 
declined  far  more 
than  the  observed 

average  HR.  Tracking  of  both  SV  and  HR  was  better  during  rapid  skin  heating  than  during 


Figure  13.  Time  course  of  selected  simulator  outputs  aixl  6  experimental  counterparts 
from  the  second  part  of  Study  No.  5. 
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rapid  skin  cooling. 


Discussion 


This  paper  documents  the  development  schemes  of  a  relatively  simple  computer 
simulation  of  thermoregulation  designed  to  predict  the  time-course  of  certain  critical  body 
temperatures  and  physiological  variables  in  response  to  the  demands  of  exerdse  and 
heat  stress.  The  capability  of  this  simulation  was  measured  by  comparing  outputs  with 
data  from  seven  independent  studies.  Overall,  the  variance  among  subjects  and  the 
variance  between  data  and  simulation  were  quite  similar,  suggesting  that,  under  the 
conditions  tested,  the  computer  simulation  may  be  able  to  forecast  an  outcome  with  the 
same  order  of  precision  as  a  laboratory  study  of  the  same  problem. 

Every  component  algorithm  in  the  simulation  is  an  approximation  and  contains  error. 
For  instance,  the  level  of  SV  is  dependent  upon  many  factors  but  only  Vp,  and  were 
included  in  the  present  algorithm;  values  for  BF,^  over  the  entire  body  were  inferred  from 
measurements  at  one  forearm  site;  values  for  are  entirely  hypothetical;  et  cetera. 
It  is  somewhat  surprising,  therefore,  that  such  good  agreement  was  obtained  between 
simulation  and  experimental  data.  Part  of  the  reason  may  lie  in  the  fact  that  errors  in 
opposite  directions  may  tend  to  cancel  within  closed-loop  systems.  In  any  ceise, 
satisfactory  overall  agreement  should  not  be  constmed  as  a  recommendation  for  the 
correctness  of  any  particular  algorithm. 

Several  differences  between  data  and  simulation  that  appear  to  result  from  systematic 
errors  in  estimation  of  regional  blood  flow  distribution  and  of  SV  during  transitions  in 
workload  and  were  pointed  out.  It  should  be  emphasized  that  the  algorithms  for 
regional  blood  flow  and  SV  were  developed  from  a  very  small  database.  These  problems 
indicate  areas  where  further  research  would  be  worthwhile. 

There  is  one  final  caveat.  This  simulation  and  its  component  algorithms  are  descriptive 
models  ratherthan  explanatory  models  (16).  That  is,  they  are  'empirical  black  boxes'  that 
are  used  to  describe  functional  relationships  between  certain  thermal  and  physiological 


responses  and  there  is  no  attempt  here  to  imply  direct  causal  relationships  between 
inputs  and  outputs. 


SUMMARY  AND  CONCLUSIONS 

1 .  This  report  documents  a  new  computer  simulation  of  physiologicai  and  biophysical 
function  designed  to  aid  in  forecasting  the  time  course  of  responses  to  various 
combinations  of  climate,  ciothing  and  exerdse  intensity.  The  body  is  modeled  as  a 
single  cylinder  containing  five  concentric  annular  tissue  compartments  (core,  muscle, 
fat,  vascular  skin  and  avascuiar  skin)  and  an  interconnecting  central  blood  compart¬ 
ment.  Algorithms,  dependent  upon  body  temperatures  and  exerdse  ievei,  control 
compartmental  blood  flows,  stroke  volume,  sweating  rate  and  shivering.  Surface  heat 
exchange  between  ciothing  and  environment  is  calculated  using  the  accepted 
techniques  of  partitional  calorimetry.  Inputs  indude  environmental  variables,  clothing 
parameters,  workload  and  anthropometric  data  Outputs  include  compartment 
temperatures  and  blood  flows,  heart  rate,  stroke  volume,  sweating  rate,  fluid  loss  and 
heat  storage. 

2.  Comparison  of  simulator  outputs  for  core  and  skin  temperatures  and  heart  rate  with 
average  data  frc'  ■  seven  independent  studies  indicates  reasonable  agreement;  on  the 
whole,  differences  between  data  means  and  simulation  outputs  were  no  greater  than 
the  data  standard  deviations. 

3.  Validation  of  the  simulation  in  general  does  not  automatically  verify  the  physiological 
control  algorithms  incoqx>rated  within  the  simulation,  since  in  complex  systems  errors 
in  system  elements  may  be  in  opposite  (Trections  and  cancel.  Systematic  discrepan¬ 
cies  in  heart  rate  prediction  and  data  in  two  of  the  validation  studies  suggest 
systematic  errors  in  either  underpredicting  cardiac  output  or  overpredicting  stroke 
volume  during  recovery  from  the  stress  of  heavy  exercise  or  of  high  skin  temperature. 

In  either  case,  it  should  be  reemphasized  that  only  a  scanty  database  based  on  a 
handful  of  studies  was  available  to  develop  the  algorithms  that  control  regional  blood 
flow  and  stroke  volume. 

I 
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RECOMMENDATIONS 


1.  Studies  are  recommended  to  determine  the  relationship  between  stroke  volume 
magnitude  with  changes  in  posture,  workload  and  recovery,  under  conditions  of  cool 
and  elevated  skin  temperature. 

2.  At  present,  the  simulation  does  not  taxe  into  account  relationships  betw-„n 
progressive  dehydration,  reduction  in  plasma  volume  and  decrements  into  stroke 
volume.  Under  certain  conditions  the  state  of  hydration  may  be  an  important  factor  in 
determining  tolerance  time,  the  level  of  cardiovascular  function  and  the  level  of  core 
temperature.  Efforts  to  incorporate  the  effects  of  hydration  status  into  the  simulation 
should  be  encouraged. 

3.  The  present  model  does  not  attempt  to  generate  a  single  index  datum  representa¬ 
tive  of  heat  stress,  heat  strain,  tolerance  or  survival  time  but  the  pattern  of  physiologi¬ 
cal  responses  accompan^ng  exposure.  It  is  up  to  qualified  professional  personnel  to 
interpret  simulator  outcomes.  'Hie  simulation  has  been  verified  by  comparison  with 
seven  human  studies.  A  larger  validation  database  is  needed.  Further,  the  present 
simulation  should  be  compared  with  other  approaches  used  to  assess  the  impact  of 
work,  clothing  and  heat  stress. 
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APPENDIX 


RECENT  ENHANCEMENT  TO  THE  SIMULATION 

A  new  feature  has  been  added  to  the  simulafion  since  the  initial  preparation  of  this 
report.  Since  this  is  an  addition  to  and  not  a  <^ange  to  the  simulation  as  described  in  the 
main  body  of  this  report,  it  is  treated  as  appendix  material. 

Hydration  Status 

The  simulation 
now  calculates 
changes  in  total 
body  water.  Ini¬ 
tially,  it  is  assu¬ 
med  that  the 
subject  is  euhydr- 
ated.  The  size  of 
the  initial  total 
body  water  com¬ 
partment  (TBW)  is 
calculated  as 
TBW  =  0.72V,*^ 

This  initial  volume 
is  simultaneously 
depleted  by  swea¬ 
ting  and  augmented  by  transport  of  water  from  the  water-load  compartment  (Rgure  At ). 
The  factor  Smiting  the  rate  of  transport  from  the  water-load  compartment  into  TBW  is  the 
rate  of  gastric  emptying,  not  the  rate  of  intestinal  absorption..  Gastric  emptying  rates  of 
15-20  mNmin'^  have  been  measured  (hiring  exercise  in  the  heat. 

An  example  of  a  recent  application  of  the  simulation  and,  in  particular,  this  new 
algorithm  was  as  an  ad  in  the  design  of  a  field  investigation  of  soldiers'  physiological 
responses  to  work  in  the  desert  at  FL  Bliss  while  wearing  protective  clothing.  Part  of  the 


Figure  A1.  illustration  of  algorithm  for  calculaling  the  time  course  of  change  in  total 
body  water. 
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protocol  required 
soldiers  to  wear 
Battle  Dress  Un¬ 
dergarments  and 
Battle  Dress  Ov¬ 
ergarments  in 
MOPP4  configu- 
rcition  while  mar¬ 
ching  at  1.12 
m-s''  with  a  22  kg 
load  for  24  min¬ 
utes  out  of  each 
30  min  period. 

Forty-five  minutes 
prior  to  the  march 
subjects  were 
preloaded  with 
450  ml  of  water. 

The  question 
arose:  "Will  sub¬ 
jects  exceed  the 
5%  limit  on  hydra¬ 
tion  before  they 
exceed  the  upper 
rectal  temperature 
limit  of  39.5°Cr' 

Two  conditions 

were  tested:  (1)  no  further  water  and  (2)  500  ml  of  wafer  every  30  min.  Panel  (a)  shows 
estimated  energy  expenditure  while  panel  (h)  shows  estimated  dry-bulb  temperature  and 
mean  radiant  temperature.  These  were  estimated  from  typicai  temparatures  in  El  Paso, 
TX  during  the  month  of  July.  The  simulation  results  indicate  that  if  the  subjects  drank  no 
more  water,  the  5%  dehydration  limit  and  the  39.5°C  Emit  on  rectal  temperature  will  occur 
at  nearly  the  same  time.  On  the  other  hatKt,  if  subjects  drink  1  ('hr'  one  would  antirapate 
hyper-hydration  to  the  extent  of  1-2%  even  when  sweating  at  a  maximal  rate. 


ri.  8tn  Ptotocoi  Time  (minutes)  scemuho  •imuiotiM 
7*  I,.  J76  cw.  sr  07/OI/Sl 

sou  »  TOO.  7«ll  -  II  AH 


Flguie  2A.  Simulation  ot  physiological  msponses  to  protocol  of  FL  Bliss  fieU  study. 
Antkapatfid  diy  bulb  and  mean  rartant  tempeiatures  were  projscled  from  other 
studies. 
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Glossary 


Af,  =  Dubois  surface  area  of  subject  or  cm^ 

=  area  of  cylinder  at  midpoint  between  adjacent  compartments  [cm] 
A,  =  effective  radiating  area  of  the  body  surface  [m^ 

BF„  =  rate  of  blood  flow  through  compartment  n  [cm^^min  '  or  t*min''] 

C  =  rate  of  convective  heat  loss  from  skin  surface  [W>m'^ 

=  rate  of  convective  heat  loss  from  respiratory  tract  [W  or 
CO  =  caidiac  output  [cm'«min'’  or  f^in'^] 

£  =  rate  of  evaporative  heat  loss  from  skin  surface  [W«m'^ 

=  maximum  rate  of  sidn  evaporative  heat  loss 
=  rate  of  evaporative  heat  loss  from  respiratory  tract  [W  or  W*m^ 
Fj  =  Burton's  thermal  efficiency  factor  for  clothing,  corrected  [MD] 

=  Nishi's  modified  permeation  effidency  foclor  for  dothing  [ND] 

H  =  upright  stature  of  subject  [cm] 

H„  =  rate  of  heat  production  in  compartment  n  [W] 

HR  E  heart  rate  [beats«min''] 

1,1,  =  thermal  insulation  provided  by  dothing  [do] 

K„„  =  conductance  between  compartments  m  and  n  [W^'C’] 

L  =  cylinder  height  [cm] 

=  total  rate  of  metabolic  energy  transformation  [W  or  W«m'*] 

PmM  =  ''apof  pressure  of  sweat  on  skin  at  temperature  [Torr] 

P„  =  ambient  water  vapor  pressure  [Tonr] 

AQ„  =  change  In  heat  content  of  compartment  n  [W>min] 

R  =  rate  of  radiant  heat  loss  from  skin  surface  [W^m"*] 

SV  =  cardiac  stroke  volume  [ml] 

T,  =  ambient  air  (dry  bulb)  temperature  [X] 

Tj  =  clothing  temperature  [X] 

=  esophageal  temperature  [X] 

T,  =  temperature  of  compartment  n  [X] 

T,  =  ambient  operative  temperature  [X] 

T,  =  mean  radiant  temperature  [X] 

T„  =  rectal  temperature  [X] 
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s  mean  weighted  skin  surface  temperature  I°C] 

V  =  subject  totai  body  voiume  [cm^ 

VqM  =  voiume  of  the  nth  annular  compartment  [cm^ 

=  subject  body  fat  volume  [cm^ 

Vrfa  =  subject  body  non-fat  volume  [cm^ 

=  volume,  cylinder  center  to  extent  of  the  nth  compartment  [cm^ 

VOj  =  oxygen  consumption  [(•min''} 

Vo^max  =  maximal  oxygen  consumption  ((•min''] 

W  =  subject  body  weight  [kg] 

W=  rate  of  energy  expenditure  dissipated  in  external  work  [W  or  W'm'^ 
c„  =  heat  capacity  (specific  heat)  of  compartment  n  [W»min^g  ’»'’C''] 
f^  =  Breckenridge  clothing  area  factor  [NO] 

h  =  combined  non-evaporat'ive  skin  heat  transfer  coefficient  [W>m'^*°C''] 
hj  =  convective  heat  transfer  coefficient  for  skin  [W^m'*.'’C''] 
h,  =  evaporative  heat  transfer  coefficient  for  don  [W^m'*»Torr'] 
h,  =  linear  radiant  heat  transfer  coefficient  for  skin  [W^m'^«°C''] 
i„  -  Woodcocl:  mcisture  index  of  water  vapor  through  clothing  [NO] 
k„  =  thermal  conductivity  of  compartment  n  (W^cm»®C''«m'®] 
rh„  =  rate  of  sweat  production  [gamin''] 
r  =  outside  cylinder  radius  [cm] 

r,„  =  radius  to  the  center  of  mass  of  the  nth  compartment  [cm] 
r„  =  inner  radius  of  the  nth  annular  compartment  [cm] 

A„p  =  area  of  cylinder  at  midpoint  between  ac^acent  compartments  [cm] 
r„^  =  radius  of  cylinder  at  midpoint  between  adjacent  compartments  [cm] 
t  =  time  [minutes]’ 

v„  =  velocity  of  air  movement  on  stationary  subject  [m^sec  ’J 
Vmo«  =  velot%  of  subject  mcvement  in  stafionary  air  [m^sec  ’] 

=  length  of  conduction  path  between  compartments  m  and  n. 
p  =  subject  average  body  density  [g-cm"^ 
p,^  =  density  of  body  fat  components  [g-cm^ 
p„  =  average  density  of  compartment  n 


'The  simulation  is  constructed  around  the  minute  time  base 
and,  unless  otherwise  specified,  values  for  variables  and  constants 
are  in  terms  of  minutes. 
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=  density  of  non-fat  body  components  [g-cm"^ 
a  -  Stephan-Boltzman  constant:  S.sy^lO"*  [W^m^^K"*) 

X  =  latent  heat  of  sweat  vaporizjtfon  at  T  •=  35®C  [W*mm*g  ’  or  J«g'’] 
%h,  =  percent  of  subject  body  weight  that  is  fat 
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